First-principles investigation of transient current of molecular devices 
by using complex absorbing potential 
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Based on the non-equilibrium Green's function (NEGF) coupled with density function theory (DFT), namely, 
NEGF-DFT quantum transport theory, we propose an efficient formalism to calculate the transient current of 
molecular devices under a step-like pulse from first principles. By combining NEGF-DFT with the complex 
absorbing potential (CAP), the computational complexity of our formalism (NEGF-DFT-CAP) is proportional 
to 0{N) where is the number of time steps in the time-dependent transient calculation. Compared with 
state-of-the-art algorithm of first principles time-dependent calculation that scales with at least A'^, this order N 
technique drastically reduces the computational burden making it possible to tackle realistic molecular devices. 
To ensure the accuracy of our method, we carry out the benchmark calculation compared with exact NEGF- 
TDDFT formalism and they agree well with each other. As an illustration, we investigate the transient current 
of molecular device AI-C3-AI from first principles. 

PACS numbers: 71.15.Mb, 72.10.-d, 85.65.+h, 73.63.-b 



INTRODUCTION 

With the advance of nanofabrication techniques, people can 
fabricate the nano-devices using single atoms or molecules 
from bottom-up approach, which leads to a new field of 
molecular electronics yj-jTi] . Many experiments have been per- 
formed to measure quantum transport properties of molecu- 
lar devices d^-lztl- At the same time, people have made a lot 
of research efforts to understand these properties from first 
principlesji-lll. At present stage, quantitative agreement 
between theoretical first principles calculations and experi- 
ment results can be reached when the system is in the steady 
state regime under external DC bias lll2l - jl5ll . Besides the DC 
steady state problem, the question of how fast a molecular 
device can turn on and off is also an important issue, which 
attracts a lot research attention recently (ll6|-|21|. This kind of 
question can be answered by studying the dynamic response 
of molecular devices by sending a step like pulse from the 
electrodes. For this problem, exact solution of transient cur- 
rent was obtained by Wingreen et al lll6ll in the wide-band 
limit using non-equilibrium Green's function (NEGF). Re- 
cently this solution has been extended to the regime of fi- 
nite band width of electrodes |{2C|. When applying this exact 
NEGF solution to molecular devices for calculation of tran- 
sient current as a function of time, there is a huge compu- 
tational cost that scales with 0{N^) due to the triple integral 
over energy, where is the number of time steps. In addition, 
there are many quasi-poles near the energy axis making the 
integration very difficult to converge. On the other hand, the 
theoretical prediction of the transient dynamics of molecular 
devices from first principles can be addressed by numerically 
solving scattering wave function or non-equilibrium Green's 
function (NEGF) combined with time dependent density func- 
tional theory (TDDFT)|[lMl- 

These methods again are very 
time-consuming for transient current calculation although the 
scaling has been reduced to N^{log2N)~ . Therefore, to speed 
up calculation, various approximate schemes were proposed 
to calculate time dependent transient current of molecular de- 



vices such as wide band approximation ^ ill . Another approx- 
imate scheme based on exact NEGF solution was also pro- 
posed and applied to calculate transient current of molecular 
devices which is very efficient and goes beyond the wideband 
limit l22ll . Despite of these efforts, time-dependent calcula- 
tion of transient current for molecular devices is still a chal- 
lenge on the computational resources. Due to the importance 
of molecular electronics, it is timely to overcome this prob- 
lem so that realistic transient dynamics calculations can be 
performed on molecular devices from first principles. In this 
paper, we propose a linear scaling 0{N) scheme to calculate 
the time dependent transient current by combining complex 
absorbing potential (CAP) method with exact solution based 
on NEGF Q and DFT theory (NEGF-DFT-CAP). 

The CAP was initially used to simulate the time-dependent 
evolution of wave function of finite systems in one and 
two dimensions [23]. Recently, CAP was employed to 
study the transport problem of molecular device from first 
principles i24l - l27ll using a transmission free CAP ll30ll . By 
adding an energy independent CAP in lead regions, the trans- 
port problem in a infinite open system can be reduced to that 
of a finite simulation region. Moreover, one can obtain an ef- 
fective wideband-like formula to calculate dc transport quan- 
tities such as the transmission coefficient. We note that in gen- 
eral ac transport properties do not assume the wideband form 
(effective self-energy does not depend on the energy) in the 
presence of CAP. Fortunately, for the step-like pulse, we are 
able to cast the exact NEGF solution for transient current into 
a wideband form using CAP which enables us to speed up the 
calculation tremendously. In fact, when CAP is implemented 
into the exact NEGF solution the amount of calculation scales 
like cNir' where is the number of time steps, « is the di- 
mension of the Hamiltonian in the whole simulation region 
including the CAP region, and c is a constant of order of a few 
hundred. With this order method at hand, the first princi- 
ples calculation of transient current of realistic molecular de- 
vices are within the reach. We have applied our formalism to 
molecular devices and carried out the benchmark calculation 
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for a one-dimensional atomic chain which agrees with the re- 
sult from exact numerical calculation. Furthermore, we have 
investigated transient dynamics of a 3D molecular device and 
calculated transient current at two different bias voltages. It 
was found that the transient current involves many time scales 
showing that the wideband limit is a bad approximation for 
molecular devices. 

The paper is organized as follows. In section 2, we will 
first introduce the formalism of complex absorbing potential 
(CAP) and briefly discuss its application in DC transport cal- 
culations. Then we will discuss how to apply the CAP to 
calculate the time dependent transient current of molecular 
device under upward step-like pulse. In section 3, bench- 
mark comparisons with NEGF-TDDFT method is presented. 
Then numerical calculation of transient current of the AI-C3- 
Al molecular device is given. Finally, section 4 serves as dis- 
cussion and conclusion part. 

THEORETICAL FORMALISM 
Complex absorbing potential (CAP) 

As shown in figure 1, a typical two terminal device consists 
of the central scattering region connected by two semi-infinite 
external leads along the transport z direction. The correspond- 
ing Hamiltonian of the whole system can be expressed as a 
tri-diagonal block matrix 





Hll 


Hlc 







H = 


HcL 


Hcc 


HcR 


(1) 







Hrc 


Hrr 





where Haa, a - L,R is the semi-infinite Hamiltonian of lead 
a. In order to study the transport properties of this open sys- 
tem, one is actually solving the scattering problem of the infi- 
nite dimension. In the framework of non-equilibrium Green's 
function (NEGF), one calculates various Green's function of 
the central region and the effect of the leads is taken into ac- 
count by the self-energy. For instance, the retarded Green's 
function of central region in energy domain is defined as 



a=L,R 



where 2J'j(£) is the self energy of lead a 

I'^iE) = Hcag'aaiEWaC, 



(2) 



(3) 



and gaaiE) is the retarded Green's function of the correspond- 
ing lead a 



(4) 



After obtaining the retarded Green's function, one can calcu- 
late various transport quantities, such as transmission coeffi- 
cient 



m 



FIG. 1: Schematic plot of a two terminal molecular device. The 
device consists of central molecular part (green solid cube) and two 
semi-infinite leads, which will extend to the ±00. The black solid 
lines represent the complex absorbing potential added to both lead 
regions. The region enclosed by red dashed line is the central region; 
the purple dashed dot line encloses the central region plus complex 
absorbing potential region in the leads. 



Here r„(£') = iCL'^iE) - I."„{E)) is the Unewidth function 
of lead a. In the numerical calculation, the energy depen- 
dent self-energy can be calculated by using the iterative or 
quadratic eigenvalue approaches 11281 12911 . To distinguish from 
the CAP method, we will refer the above method as exact 
method. 

The idea of CAP method is to replace the infinite system 
by a finite dimension using a transparent boundary condition 
that absorbs the incident wave function completely. In the ap- 
plication of CAP to the quantum transport problem, CAP is 
added to a finite lead region (called CAP region) outside of 
central scattering region. Usually, the effectiveness of CAP 
on absorbing the incident wave depends on the length of the 
CAP region. The reduction of reflection can be improved by 
increasing the length of CAP region in a controlled way. Note 
that the advantage of CAP method over the exact method re- 
lies on the fact that the CAP does not depends on energy while 
the self-energy of the exact method does. Using this property, 
the poles of Green's function for CAP method can be obtained 
easily. Therefore the convergence problem of energy integral 
in calculating transient current is solved. In the numerical cal- 
culation, we adopt a typ e of optimized transmission-free CAP 
form given in Ref. OOll . 



Wiz) = —i—ffiz), 
2m Az 



where f{z) is defined as 



4 Az 

f(z) = -(( 

c- zi -2zi + z 



f + i^f-2), 
zi-z 



(6) 



(7) 



T{E) = Tr[riG^cr«G^c]- 



(5) 



and Az - Zi - Z\ is the range of CAP along transport z di- 
rection, z\ and Z2 are the starting and ending points of CAP 
region at each lead, respectively. Here c is a constant taken 
to be 2.62, m is the electron's mass. As shown in figure [1] 
the CAP region starts from several buffer layers away from 
the central molecular region. Going deep into the lead, the 
strength of absorbing potential increases and /(z) 00 when 
z approaches the end point zi- This truncates the semi-infinite 
lead into a finite one. Therefore, the numerical simulation re- 
gion becomes finite, i.e., the region enclosed by the purple 
dashed dot line in figure 1 . In the first principles calculation. 
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LCAO basis set is usually adopted. Then one has to calculate 
the matrix element of CAP in orbital space 



y, z)W{z)(f>v(x, y, z)dxdydz. 



(8) 



where 0^ is the atomic orbital. 



Within the CAP method (all quantities are labeled with 
prime), the retarded Green's function of lead a can be defined 
asS 



(9) 



Since the lead is effectively truncated, H'^^ - Haa - i^a is 
a matrix of finite dimension shown in figure 1 . The retarded 
Green's function of whole system including the CAP region 
can be expressed as 



C' + / J] C)" 



(10) 































and W'^ = 





























with W[ 



Although the lead region (CAP region) is finite, we can still 
use the concept of self-energy and obtain an eff'ective retarded 
Green's function of the central region within CAP method 



'cc 



(E - Hcc 



where self energy EJ^ is given byyj. | 

i;;'(£) = HcaS'LiEWaC- 



(11) 



(12) 



It is easy to show that the linewidth function is written as 113211 

fff - 2Hcag'aa^a8aafiaC = '^Hcagaa^ag'aaf^aC ■ (13) 

Since the self energy calculated by CAP method is the same as 
that obtained from the traditional method[26], various Green's 
functions in the central (physical) region shown in figure [T] 
should also be the same as that given in Eq. (|2]i. At this 
level, the self-energy of the Green's function of the central 
scattering region GJ^^ depends on energy. In the following, 
we give a simple derivation on transmission coefficient us- 
ing CAP method in an effective wide band limit (WBL) form. 
Starting from the traditional definition of transmission coeffi- 
cient of Eq. (fT2l i together with Eq. (fT3] l 

TiE) = Tr[rlG^cr«G"cc] 

= 4TT[HcLgi,WLgiLHLcG'ccHcRgiRWRgi^HRcG''cc] 
'fJl^^TrlW'^C'W'ifG"'] (14) 



4Tr[WLG',^WuGl^, 



where we have defined the following Green's function of the 
whole system including the CAP region (see the Appendix A 
for derivation) 



In order to calculate transmission coefficient of Eq. (fl4] i. one 
only needs to know which is defined in the whole sys- 
tem including the CAP region with W', an effective energy 
independent self energy. Note that this effective WBL form is 
only valid in DC case. In the case of AC transport, one may 
not have similar WBL form and one has to deal with it case 
by case. 

In terms of the lesser Green's function, one can calculate 
the charge density in the central region. We also start from the 
traditional definition of the lesser Green's function 

Gcc(E) = i ^ faGccKG'cc 

a 

= 2iY,UG'c,'^aGic 

a 

^2iJ]fa[G'"w:fi'']cc, (16) 

a 

where we have used Eq. (fT3T l and GJl.^ = G''^fjHcagaa (sse 
Eq.(3.5.13) in Ref. i3lll ): fa is the Fermi distribution function 
of lead a. 



Time dependent transient current with upward step like pulse 

The exact solution of time dependent current for step like 
pulse based on NEGF has been given by Maciejko et al ll20ll . 
This formalism can be combined with DFT to calculate tran- 



sient current in molecular devices Il22[ 13311 . In the following, 
we will combine the exact solution with DFT and CAP to ob- 
tain an order 0(N) scheme (NEGF-DFT-CAP) for calculating 
time dependent current under upward step like pulse. Down- 
ward step and square like pulses can also be treated in a similar 
fashion. 

To begin with, we will derive an equivalent time depen- 
dent current formula. Startjrig from the equation of motion 
for lesser Green's function 113411 . we have 



Hcc{t)G^cc(t't')+ f [^^(t,ti)G''cc(tut') 
Jo 

+I.''{t,ti)G''{ti,t')]dti, (17) 



and 



G'lr - g'LL^LcG'ccHcRgRR- 



(15) 



-i^G^c(t't') = G^c(t't')Hcc(t') + j [G[.cit,h)^^{tut') 

+G''{t,H)If{tut')]dh. (18) 

Then subtracting equation (fTST l by equation ( fTTj i and setting 
t' - t, we can arrive at 

lop(t) = Hcc(t)G^c(t' t) - Glcit, t)Hcc(t) - ijG^cit^ tU19) 
where we have defined 

lopit) = r [G'cc(t,ti)^-'(tut) + G^{t,ti)I."{ti,t) 
Jo 

-S<(f, ti)GlcitiJ) - S'Xf, ti)G\ti,t)]dti. (20) 
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which is a matrix. Note that the terminal current /n.(f) l35ll can 
be obtained from Iop{f)- To do that, two auxiHary projection 
matrices are introduced 



(21) 



where a denote outermost unit cell layer in the buffer layers 
of central region and is the unit matrix with dimension 
equal to the size of unit cell of left and right lead respectively 
so that = f ^Z'^f i. Finally we have 
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h(t) = Tr[r„/„„(f)r„] 



(22) 



From ( fT9] l, we can use the following formula to calculate the 
time-dependent terminal current /„(f), 

hit) = 2R&Tx[Y,Hcc{t)Glcit,tta] - /Tr[r„5,G^c(f, Ofj, 

(23) 

where G^f.{t, t) is the time dependent lesser Green's function 
of central region with equal time. To calculate lait), one has 
to know the time dependent Hamiltonian H{t) and calculate 
time dependent lesser Green's function G^f.{t, t). Since the 
external bias is the upward step like pulse in our problem, then 
the time dependent Hamiltonian can be obtained as follows. 
When time f < 0, H{t < 0) = Heq is equilibrium Hamiltonian 
without bias and H(t > 0) = 0{t)H„eq that H„eg is the self 
consistent non-equilibrium Hamiltonian under DC bias. As 
for the time dependent lesser Green's function within CAP 
method, it can written as 

G<c(t,t)^2iJ] J ^f(oj)[A'^(oj,t)W'J^{oj,t)]cc, (24) 



where we have used the spectral function A^{cij, m35\ 

A'„(ej)= f t/f'e''^"-''>e'£'"'^"("'G'''(f,f'), (25) 

«7— oo 

where A„(f) is the time dependent external bias. Note that 
A^(e, f) has the same dimension with G' that is defined in 
CAP method. Then the key issue is how to calculate quantity 
A'g,{e,t) efficiently. From the analytic expression of A„(e, f) 
given in Ref. 12011 . we can derive the spectral function in the 
CAP form (see Appendix for derivation), 

ddj e-'("-^>'G''(w + A„) 



A>,f) = G''(e + A,)- J 
A 



■iO 



2ni cl) — £ + Aa — iO^ 
7+AG^'(e)], 



(26) 



where A - Hn^q - Heq is internal potential change due to 
the external bias. Equilibrium and non-equilibrium retarded 
Green's function are defined as 



G''{e)^[eI-H,q + iY^W'^r\ 



(27) 



FIG. 2: Schematic diagram of a molecular device AI-C3-AI. The 
device consists of three carbon atoms chain coupled to the perfect 
aluminium atomic electrodes which will extend to the reservoirs at 
±00, where the current is collected. 



G''(e)^{eI-H„,q + iY^W'S'- 



(28) 



Since is energy independent, we can use following eigen- 
equations to construct the retarded Green's functions, 

[H,„,-/^W;]|.A«) = e»l«An) (29) 

a 

[Hneq + /^(C)"'"]l0n> = ^l-^"). (30) 

a 

and 

[H,q-iY^W'M)^(^Mr) (31) 

a 

[Heq + i YJ-W'JmI) = e°*|^°>. (32) 

a 

Then retarded Green's functions can be constructed from their 
eigenfunctions 



{e-€i; + /0+) 



)■ 



(33) 



(34) 



Due to the presence of the time dependent factor g 'f*^ f)' 
in A„, the integration in A^(e, t) can be done analytically by 
enclosing a contour in the lower half of complex plane. 



A>,o=y — — +y 

" 4^ (e + A„ - e„ + ;0+) ^ e-e„ + /0+ 



/(e+A„-e„)/ 



\^n){4>n\ 



e + A„-e„ + /0+ Ai(e_gO+/o+) 



(35) 



It is easy to check that in the initial and asymptotic long 
time limit (f — > oo) A^(f) is equal to G' (e) and G' (e), re- 
spectively. After obtaining the A^j(0, one can calculate lesser 
Green's function using Eq. (|24] | and hence the time depen- 
dent current from Eq. (l2Jt . Since A^(f) is expressed as the 
summation form at any given time t, one only needs to inte- 
grate the energy co in Eq. (l24l i to obtain the time dependent 
lesser Green's function and hence the transient current /„(?)• 
We can estimate the number of operations in calculating time 
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FIG. 3: The time dependent transient current I(t) versus time with 
V = 0.0005 a.u. for one dimensional atomic Al-Ci-Al chain. Blue 
solid line and red dashed line are time transient dependent current 
calculated by using NEGF-TDDFT method and NEGF-DFT-CAP 
method; black solid line is the DC current at steady state limit. 



dependent current /„(?)■ For a given time, the calculation only 
involves matrix multiplication as well as the integral over lj 
which again can be done using the theorem of residue. Hence 
the total number of operations is roughly cNn^, an order N 
algorithm, where c is of order 200 due to the contour integral 
on the complex plane, « is the dimension of Green's function 
of the whole system including CAP region, and is number 
of time steps. 

The major steps for the numerical calculation can be sum- 
marized as follows. We first prepare the initial equilibrium 
and final non-equilibrium self consistent Hamiltonian from a 
DC calculation. Then we construct the CAP matrix Wa with 
respective to the lead. Once the CAP is constructed, one can 
compare the transmission coefficients with that obtained by 
exact method to get an idea how long the CAP region should 
be. With the good agreement on the transmission coefficient, 
we can move on to calculate the time dependent current using 
A'M 

It is worth mentioning that, in the above discussion, the or- 
thogonal basis set is implicitly assumed to expand the Hamil- 
tonian. So one has to orthogonalize the basis set if non- 
orthogonal basis such as atomic orbital basis set (LCAO) is 
used.il 



NUMERICAL RESULTS 

In this section, the implementation of our formalism and 
numerical results of transient current for AI-C3-AI molecular 
device will be presented. The structure of AI-C3-AI molecular 
device is shown in figured There are 75 atoms in the central 
scattering region and the distance between the Al atom and the 
nearest carbon atom is equal to 3.78 a.u. As for the electrodes, 
there are 9 aluminum atoms in a unit cell with a finite cross 
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FIG. 4: The comparison of transmission coefficient of a Carbon chain 
sandwiched between Al(lOO) leads. The numerical results calculated 
by using CAP method with 30 unit cells in the lead region (blue solid 
line) is compared with exact method (red dashed line). 



section along (100) direction in the semi-infinite aluminum 
lead. 

Our numerical analysis is based on the state-of-the-art first 
principles quantum transport package MATDCAL. 11361 |37 I 
Specifically, a linear combination of atomic orbitals (LCAO) 
is employed to solve KS equations. The exchange-correlation 
is treated at the LDA level and the nonlocal norm-conserving 
pseudopotential l38ll is used to define the atomic core. The 
density matrix is constructed in orbital space and the effec- 
tive potential is obtained in real space by solving the Poisson 
equation. The accuracy in the self-consistent iteration is nu- 
merically converged to 10""^ eV. The initial equilibrium and 
final non-equilibrium Hamiltonians are prepared using MAT- 
DCAL package. 

In the following, the case of upward step-like pulse (Viit) = 
-Vuit) - 0{t)V) applied on both leads will be considered. In 
order to satisfy the current conservation condition, we will 
plot the time dependent current in terms of /(f) - [Iiit) - 
lR{t)]/2M- 

Before presenting our numerical results, we have calculated 
the transient current using two different approaches to test the 
accuracy of our present scheme. One is based on the NEGF- 
TDDFT method proposed in Ref. llSll which is an order 
N^{log2N)- algorithm and other one is our proposed formal- 
ism in this paper termed as NEGF-DFT-CAP. Here we take 
one dimensional Al-Ci-Al atomic chain (where both leads are 
one-dimensional Al chain) as a toy molecular device and ap- 
ply a step-like pulse to test numerical implementation of our 
formalism. As shown in figure |3] transient current calculated 
from two different methods agree well with each other In 
addition, the transient current approaches the DC steady state 
value obtained by using Landauer-Btittiker formula in the long 
time limit. The insert figure shows the early time behavior of 
the transient current. 

Now let us study the real molecular device AI-C3-AI. First 
of all, we have to compare the transmission coefficient by us- 
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Derivations for G', 



FIG. 5: The time dependent transient current /(/) versus time with 
different bias voltages for AI-C3-AI molecular device. The blue and 
black solid lines correspond with time dependent transient current 
and DC current at steady state for V = 0.0025 a.u., respectively. The 
red solid and black dashed lines are time dependent transient current 
and DC cuiTent at steady state for V = 0.01 a.u., respectively. 



According to the definition of retarded Green's function. 
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ing CAP method and exact method to make sure that the CAP 
potential is added correctly. As you can see in figure |4] the 
CAP results agrees well with that calculated by exact method. 

After comparing the accuracy of transmission coefficient of 
CAP method, we are ready to study the transient current of 
AI-C3-AI device. We calculate the transient current under two 
different bias voltages. The numerical results are plotted in 
figure |5] We have serval observations: (1) the switch-on time 
is roughly 2 fs; (2) the relaxation time is roughly 210 fs for 
V = 0.0025a.u. and 320 fs for V = O.Ola.u.; (3) the tran- 
sient current is on the same order of magnitude as that of the 
DC steady state limit. In the early time, there are some ir- 
regular oscillations in the transient current. At the long time 
limit, the transient current approaches to the correct DC limit. 
Moreover, more oscillations occur with the increase of bias 
voltage. The nature of the oscillation can be attributed to the 
resonant states of the svstem ll22ll . 



we have 



G'lR - glL^LcGcR- 



(36) 



(37) 



To find G[^^ we note that the advanced Green's function can 
be obtained by replacing superscript r into a in equation ( [36l ). 
We have 



G^RC - 8rrHrcG"cc' 
G'cR = G^cHcRgRR, 



(38) 



where we have used the fact that GJ^^ = (G^p)^ Finally we 
combine equations ( l36b and ( |38] ) to arrive at 



G'lr - glL^LcGccHcRgRR- 



Derivations for A' 



(39) 



CONCLUSION 

To summarize, we have proposed an order N first-principle 
formalism to study dynamical response of molecular device 
due to the time dependent step like external bias. Our for- 
malism is based on NEGF combined with DPT as well CAP 
method. The use of CAP allows us to calculate transient cur- 
rent efficiently. Comparing with the previous NEGF-TDDFT 
schemes, the computational cost of current scheme is much 
less. Detailed computational procedures for first principles 
transient current calculation were discussed which is very 
easy to implement. As an illustration, we have calculated 
the transient current of AI-C3-AI molecular devices from first 
principles. 



In this appendix, we will derive the expression of A'^^^ in 
central region within CAP method to confirm A[-^ given in 
equation ( [26] l. Starting from equation (|26] l, we have 



d(jj 



u - e - /0+ 



27r/ (jj - e + Aa — /O^ 
BB], 



(40) 



with 



BB= Yj G'cpioJ + K)AppG'l,c{e) 



I3=L,C,R 



(41) 



where Acc is the internal potential change in the central region 
due to the external bias in the leads and A^^ = ^p^p with 
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j3 - L,R is the bias applied in lead /?. Furthermore, BB can be 
separated into two parts. 



BB - G'(j(.AccG'(^c + ^P^cfi'pc 



P=L.R 



P=L,R 



(42) 



cc 



where we have used G[,p 



G'HcBSoa and GL 



Spp^pcG'cc- According to the retarded Green's function of 
lead a in equation (|9]l, we have 



1 



1 



1 



(43) 



Therefore, becomes 



BB - G'ccAccGcc ~ ^/^^'cc^a/si^- (^)Gq 

P=L.R 



where we have defined 



£J(6)-£J(w + A„-A^) 
e — (jj - Ac, + Ar 



(44) 



(45) 



Finally, by plugging equation (l44l i for BZ? into equation ( l40b . 
one can easily find that the final expression for A^,^^ is the 
same as the expression for Aa^c given in Ref. ll20ll . 
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